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Abstract 

A nonlinear coupled system descriptive of multi-ion electrodiffusion is investigated and all param- 
eters for which the system admits a single-valued general solution are isolated. This is achieved via 
a method initiated by Painleve with the application of a test due to Kowalevski and Gambler. The 
solutions can be obtained explicitly in terms of Painleve transcendents or elliptic functions. 

PACS: 02.20. Qs, ll.lO.Lm 

1 Introduction 



J^ ' The theory of electrodiffusion has its origin in the liquid-junction theory of Nernst and Planck [1]. It 

H I seeks to provide a macroscopic description of the transport of charged particles through material barriers. 

AppHcations abound, in particular, in the modelling of biological membranes [2]- [7] and in electrochem- 
istry [8] . Schlogl [9] observed that it is convenient to partition the ions into m classes characterized by the 
same electric charge qj = q^Vj, where go is the unit of charge and Vj is a nonzero integral signed valency. 
The TO-ion electrodiffusion model in steady regimes then reduces to the following system of tti -I- 1 coupled 
first order ODEs [10]: 



— — = ViUiP ~ Ci, ViTli ^t 0, Z = 1, • • ■ , TO, 

ax 
dp 
dx 



y^^Vjni, 



(1) 



i=i 



where x is the coordinate normal to the planar boundaries, p is the electric field and rij is the number of 
ions with the same charge qj = q^Vj. These variables are subject to the constraint [9] 

{vi - Vj)vini 7^ 0, i ^ j. (2) 
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The system (1) admits the first integral 



K 



'-,-o^-Y. 



Ui, 



(3) 



where K may be set to zero without loss of generality whenever C — ^ c^ 7^ 0. Note that certain 

boundary value problems indeed require that K be zero. Here, this constraint is not imposed and it 
is demonstrated in an algorithmic manner that the electrodiffusion model for tti < 4 admits solutions 
given in terms of Painleve transcendents and elliptic functions. The Painleve reduction obtained in [11] 
is retrieved as a particular case. 

2 Single-valued solutions. A classical procedure 

For nonlinear systems such as (1), dependent on certain parameters J^i, q, there is a classical method [12] 
for determining all those parameters for which the general solution of the nonlinear system is single- valued 
around its movable singularities. An alternative approach due to Kowalevski [13] and Gambler [14] is 
also available, and all the details which are sometimes missing in these classical authors can be found in 
[15] or in summer school lecture notes [16]. 

The main steps, to be performed in the sequel, are as follows: 

1. A scaled version of the nonlinear system (1) is introduced for new fields Ni{X),P{X), via the 
transformation (this is a particular a-transformation [12]), 



X — eX, p — e ^P, rii — Ni. 



(4) 



Thus the limit e ^ generates 



dX 



v,N,P, 



V^^Q, 



l,---,m. 



dP " 

i=l 
p2 ™ 



AX 







(5) 



This canonical system results from (1) and (3) by setting q = 0, i^ = and is termed the simplified 
system, as opposed to the original com,plete system (1). If the complete system has a single- valued 
general solution, it is necessary that the simplified system has a single- valued general solution. 

2. Necessary conditions on the parameters i>i for the general solution of the simplified system to be 
single- valued are determined. These conditions are shown to be sufficient for the problem at hand. 

3. The necessary conditions on the parameters Vi^Ci for the general solution of the complete system 
to be single-valued are obtained and the corresponding PII and other reductions listed. 



3 Constraints on the charges Vi 

The transformation (4) implies that P{X) admits movable simple poles, 

P = X-' [Po + 0{x)] , x = X-Xo, Poy^O. 
A straightforward integration of (5)i reveals the leading order behaviour 

f P = x-'[Po + 0{x)], 

p2 ™ 

= -j-x-' [1 + Oix)] Y. x'^-'^" [^^ + <^ix)] 



(6) 



(7) 



Here, the Ai are undetermined constants. Since the Vi are assumed distinct (see (2)), the term in x^^ in 
(7)3 can only be matched by exactly one of the terms involving x'^^^" ■ Denoting j this matching index, 
this yields the desired value of Pq, 

Po = -^, Aj^X. (8) 

In view of the possible vanishing of the remaining constants Ai, the necessary condition that each Ni{X) 
be single-valued near X = Xq is not easy to enforce and hence we proceed directly with the other 
conditions. 

An additional necessary condition is that the linearized version of the simplified system near the local 
behaviour (7) also has the property that its general solution is single- valued near X = Xq. Since this 
linearized system is of Fuchsian type [17, Chap. XVI] near X = Xq, it is necessary that all its Fuchs 
indices be integer. This condition is here readily enforced with knowledge of Pq. Thus, let us introduce 
the function W by 

W 
P = -2— — . (9) 

By construction, W has a simple zero near X = Xq (it can be chosen such that 

limx^Xo W/{X — Xq) = 1). The integration of the simplified system then reduces to the integration of 

a single ODE for W{X), namely 



W'^^^Yl ^^ W^'"'"'/"^ = 0, (10) 

where 



2 



N.^kiW-^"'/"', (11) 

and the ki are arbitrary constants. Since limx^Xo W/{X — Xq) — 1, Laurent series in X — Xq are also 
Laurent series in W. The values of the Fuchs indices arc directly read from (10), and arc given by 

Fuchs indices = — 1, 2 — 2—, i — 1, ... ,111. (12) 



Finally, since j can be arbitrarily chosen, the diophantine condition to be solved is 

The latter condition is restrictive and, using the shorthand notation vi : i>2 '. 1^3 — a : b : c to indicate 



Vj Vi : -2— = integer. (13) 



— = —- = -^ — arbitrary, (14) 

a b c 



its only distinct solutions are 



m = 1 : no restriction, 

m = 2 : j/i : j/2 = 1 : -2, 1 : -1, 1:2, 

m = 3 : i^i : 1^2 : 1^3 = I : -2 : -1, 1 : -2 : 2, (15) 

m — A : vi : 112 '■ V3 '■ Vi = t : —2 : —1 : 2, 

?Ti > 4 : no solution. 

The necessary condition (13) is also proven to be sufficient. Thus, in the first order ODE (10) for 
W{X), all the powers of W are integers between and 4 for the values (15) so that this ODE belongs 
to the binomial type studied by Briot and Bouquet [18, pages 58-59] and its general solution W{X) is 
single- valued. The explicit expressions for P and Ni in terms of W prove the single-valuedness of the 
general solution of the simplified system (5). This completes the second step as set down in Section 2. 

Let us remark that the Fuchs indices (12) comprise a subset of {—2, —1, 1, 3, 4, 6}. 



4 Constraints on the q 

The third step involves a well known test due to Kowalcvski and Ganibicr. In view of the numerous 
possible dominant behaviours (11) of the n,, it proves convenient to eliminate the m variables Ui to 
obtain the 77i-th order ODE obeyed by the electric field p{x). Thus, taking the derivative of (1)2 to — 1 
times, one generates a van der Monde system for the m variables n^, namely [10] 



E-^ 



-^Q) 



i," 



ai^p , paq = aq_^+'^vl ^Ck, q^2,---,m. 



(16) 



fe=i 



In view of the constraint (2), the solution rii of this system is unique. Thus, the question of single- 
valuedness of the general solution {ni,p) of the complete system (1) reduces to that of p alone. The ?Ti-th 
order first-degree ODE for p{x) is defined by the determinant 



1 



1 p72 - {Cx + K) 



Ol 



= 0. 



(17) 



For the admissible values of to, as set down in (15), one obtains [10]: 



m = 


= 1 : 


p' + Vl 


\ P' 1 


= 0, 


m - 


= 2 : 


~ J 

p" - {Vl +l'2)PP +^1^2 


rp3 1 

^ - (Cx + K)p 


m - 


= 3: 


PP"'-P 


'p" - (Ui +U2 + 


^z)p^p" + {^2Vz + 



with the notation 



+ V1C1 + V2C2 = 0, 



„3„f 



+ viV2V^ -?- + {Cx + K)p^ 

- [{1^2 + l^a^lCl + (VS + Vl)v2C2 + [vi + V2)v3C3]p'^ 

- {civi + C2V2 + czV'i)p' == 0, 



(18) 



p'p"" - Zpp'p'" + -ip'p" 



PP 



+ s,{-pY" + P'p'p") + S2pV' - sa/y + '-^p' 
+ A^p^ + Ai(3p'' - pp") + A2pV - 0, 



Si{Cx + K)p^ 



Si = Vi + V2 + J^3 + i^4, 

52 = 1^11^2 + l^lVZ + l^l^A + V2VZ + 1^2VA + V?,Vi, 

53 = V2VzVa + VlVz^A + VxV2V^ + VxV2Vj,^ 

54 = VxV2V'iV^, 

M == (C2 + C3 + C^V2V'iV^ + (Cl + C3 + C^VxVzV^ 
+ (Cl + C2 + C^VxV2V^ + (Cl + C2 + 03)1^11^2^^3, 

M = (Cl + C2)vyV2 + (Cl + Cz)v\Vz + (ci + c^)vxv^ 
+ (C2 + C3)l^2J^3 + (C2 + C4)j^2i^4 + (c3 + C4)i^3'^4, 

Ai = CxVx + C2l^2 + C3i^3 + C4Z^4. 



(19) 



The next set of necessary conditions arises when one investigates the existence of the m Laurent series 



(one for each chosen j) whose first term is (sec (6)) 



+00 



^Pkix - Xq) 



fc-1 



Po 



fc=0 



2 



J = 1, 



(20) 



The locations k at which arbitrary constants may appear in the series (20) are identical to the Fuchs 
indices listed in (12). 

Whenever k reaches one of the positive Fuchs indices in the list (12), to eliminate the possibility 
of a movable logarithm in the expansion (20) destructive of the Painleve property, conditions must be 
imposed on the q, i — 1, • • • , tti. The results of this classical computation are (Appendix): 



m — 1 : 
m = 2, 
m ~ 2, 
m — 2, 
m = 3, 
m = 3, 
m ~ A, 






i^2 = 1 
V2 = I 
V2 = I 
V2 ■ V3 
V2 ■ V3 



-2 : 
-1 : 
2 : 

= 1 : - 
= 1 : - 



1^2 '■ 1^3 '■ 1^4 = i ■ —2 : —1 : 2 



no additional condition, 

Cl =0, C2 = 0, 
no additional condition, 
no additional condition, 
ci = 0, C2(2c2+3c3)-0, 

Cl =0, C2 = 0, 

Cl — C2 — C^ — C4 — 0. 



(21) 



The value K of the first integral is unconstrained. 

Whenever —2 is a Fuchs index of some family, movable logarithms may occur as well from this index 
—2, but the classical computation recalled in the Appendix cannot detect them. In order to perform this 
detection, an analysis up to the perturbation order n — 6 has been conducted in the manner described 
in [19]. All such additional necessary conditions turn out to be identically satisfied. 

In order to complete the third step, it remains to examine whether the general solution {ni,p) of the 
complete system is indeed single- valued for the cases set down in (21). It proves convenient to make use 
of one of the m variables w involved in the singular part transformation (see (9)), 



P 



-2- 



VjW 



(22) 



namely the one where \vj\ is the greatest of the tti charges. To each vanishing Ci, if any, there corresponds 
a first integral 



riiW 



2vi/vj 



^0, 



(23) 



where the exponent 2vilvj is an integer. For each pair {ci,Ck) with c^ = c^ = one can eliminate w 
between the two associated first integrals to obtain a first integral rational in [11,1,71^), namely 



2uj/vi — 2i/j/j/fe 



= constant. 



(24) 



5 The m-ion cases: m = 1, • • • , 4 
5.1 m = 1 

In this single ion case. 



ni 



\-{Cx + K), 



where w satisfies the Airy- type equation [20], 



2w' 



(25) 



w" - -^{Cx + K)w = 0. 



(26) 



5.2 m = 2 

Here, the rii, i — 1,2 are given by 

p' -^i{\- Cx -k) p'-yJl^-Cx-K 
ni = ^ ^, n, = ^ ^, (27) 

and the second order ODE for "pix) has degree one in -p" . All such second order ordinary differential 
equations which have the Painlevc property have been classified by Painlevc [12] and Gambler [21]. 
The first case 

vx:v2 = \: -2, ci ==0, c^^ 0, / + i^ipp' - vlp^ + 2Kvlv = (28) 

belongs to Class 10 of Gambler and its general solution is elliptic with 

2 

if = 0. (29) 



v = 


2" , 

V2W 


111 = 


= kiw, 


712 = 


-k2W-\ 


W 


- kiw - 


- k2W 


2i^fw^ 


The second case 



















i/i : 1/2 = 1 : -1, p" - ^lyfp^ + vliCx + K) + v^{cx - C2) = (30) 



may be reduced to the Painleve II equation 



i2 rr 

— ^ = 2C/3 + XC/ + a (31) 

under the scaling transformation 

p=— t/, 2;== — -—, k^ = -v^C, 2a = — - — . (32) 

The third case 

vx:v2^\:2, p" - Svipp' + lyfp^ - 2vf{Cx + K) + iyi{ci + 2c2) = (33) 

belongs to Class 5 of Gambler and is linearizable with 

2w' 
p= , w'" ~2iy^(Cx + K)w' -ui(ci+2c2)w:=^0. (34) 

1^2W 

5.3 m = 3 

The three cases to be considered are: 

(3a): i^i:u2:u3 = l:-2:-l, (ci, C2, cg) - (0, 0, 1)C, ' 

(3b): 1^1 -.1^2 ■■'^3 = I- -2: -I, (ci, C2, cg) - (0, 3, -2)C, I (35) 

(3c): 1.1:^^2:^^3-1: -2:2, (ci, C2, cg) = (0, 0, 1)C. , 

The class of equations to which the ODE (18)3 belongs, namely 

—p p" + bpp'p" + dp' + hp' p" + kp p + qp p + fp (36) 

+ subdominant terms = 0, 

has been investigated in [22]- [24], but the results therein are insufficient to cover the above three cases. 
For C — K — Q, the integration by elliptic functions (see (10)) was known to Chazy [22]. For arbitrary 
values of C and K^ we proceed directly with the original system (1). 



In the first case, two q vanish, so that p obeys a second order ODE and a birational transformation 
exists between p and w'^ with 



(3a): < 



p = —2 , ni = kiw, 712 — k2W , 

J' ^'^ , ,.2 f OL ,.,2 /^- , r^N- , ^2 



«;"--— + z/f -2fcii// - (Cx + K)w + — = 0, 
2w \ w J 

2 -2 p'^ - t^iV/2 + i^i^CCa; + K)p - v^C 

ik2V^W = . 

P 



(37) 



The ODE for w belongs to Class 34 of Gambier and its general solution is an algebraic transform of the 
PII function t/ (31) with 



(3a) 



o dt/ , , :/? , X K 

kivlw^ek-- + k^U^ -^(Cx + K), x = — -—, 
ax 2, k C 



-.?C, Sklk2^ ^'^+f'^\ e- 



(38) 



1. 



In the second case, since only one of the Ci vanishes, the ODE for w is of third order and a birational 
transformation links p and w with 



-2 , rii = kiw, 

l'2W 



(3b): <^ 



P 

w'" + vl [-Qkiww' - 2{Cx + K)w' - ACw] = 0, 

QkiViW = 



(39) 



2 _p" + Si^ipp' + vIp^ ~ 2vl{Cx + K)p - AviC 



P 
The ODE for w belongs to Class XIII of Chazy [22] and admits the first integral 

/2 



K2^{w + [Cx + K)k^^) w" 



(3b): <^ 



Cw' 



w 



2kiviw^ - Avi{Cx + K)w^ - 2vf{Cx + Ky — , 



ki 



(40) 



= vik^ ni [i'i{Cx + K)n3 — Cp — Vinin2] 



This second order ODE for w again belongs to Class 34 of Gambier and it may be integrated in terms of 
the PII function U with 



(3b): 



kiviw = ek- h/c^C/^, 

da; 

x^^-^, k^ = 2ufC, kJK2 - ~2C^a{a + e), e^ = 1. 
k C 



(41) 



The third case, for which again two Ci vanish, follows the same pattern as the first case and we have 



(3c) 



W 7—1 J 2 

p ~ —2 , ni — kiw , n2 — k2W , 

w" + vl [~Ak2W^ - 2{Cx + K)w - fci] == 0, 



(42) 



I2k2vlw 



2 _p" - Svipp' + i^tP - '2iy{{Cx + K)p + 2viC 



P 



so that a PII equation is obtained for w. 

It is remarked that, in the three cases considered above, the first integral of the third order ODE 
(18)3 for p is a polynomial of third degree in p" . This underlines the advantage in proceeding with w 
rather than p in our analysis. 



5.4 m = 4 

In the unique case to be integrated, 

j/i : z^2 : t^3 : J^4 = 1 : -2 : -1 : 2, Ci = C2 = C3 = C4 = 0, 
the four first integrals immediately yield the solution, 



(43) 



P = 



2l^2W 



rii — kiw, n2 = k2W , n^ = k^w , n^ — k^w 



(44) 



2vlw^ 



— kiw — k2W ^ — k^w ^ — kiuP' — K = Q^ 



so that w in general is an elliptic function. 

The case when m is arbitrary and all the c^'s vanish is similarly reduced to a first order ODE intcgrable 
by a quadrature, however with a multi-valued solution when to > 4. 

6 Conclusion 

The analysis presented here has been used to isolate, in particular, underlying Painleve II structure in 
the 3-ion electrodiffusion model. This has potential application to the construction of upper and lower 
solutions to two-point boundary value problems [25] . The intcgrable nature of the model in these cases 
also allows the application of Backlund transformations [26]. 
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Appendix. Conditions for the absence of movable logarithms 

After the coefficient po of the series (20) has been computed, the recurrence relation for subsequent p^, 
namely 

Vfc > 1 : Ek= P{k)pk + Qk{{pi \l<k})^0, (45) 

depends linearly on pk and nonlinearly on the previously computed coefficients pi . Whenever the positive 
integer fc is a Fuchs index r of the TO,-th order ODE (17), the coefficient P{r) vanishes and the condition 
Qr = must be enforced in order to avoid movable logarithms. For the cases isolated in (15) and for 
each of the m families po — ^2/^'^, these conditions Qr ~ deliver the following. 



m 



2: 



vi : iy2 = I 



vi : V2 
vi : V2 
vi : V2 
vi : V2 
vi : U2 



-2, i = 1, (36 = C2(ci 

-2, J = 2, Q3 = ci, 
-1, ] = 1, Q4 = 0, 
-1, J = 2, g4 = 0, 
2, j = 1, no condition, 
2, ) - 2, Qi - 0. 



2C2), 



(46) 



3: 



vi : V2 : V7i = '^ 

vi : 1^2 '■ J^3 = 1 

1^1 : 1^2 : vz = 1 

j/i : z^2 : J^3 = 1 

j/i : z^2 : v^ = 1 

vi : V2 ■■ V3, ^ I 



-1, i = 1, (34 = 0, Qe = C2(ci + 2c2 + 3C3) 

-1, j = 2, Qi=0, Q3=ci, 

-1, j = 3, Q4 = 0, 

2, j = 1, Qe = C2(5ci + 10c2 + 6C3), 

2, j = 2, Q3=ci, g4 = 0, 

2, j = 3, Qi = 0, Qa = C2P1. 



(47) 



m = 4: i>i : 1^2 '■ 1^3 '■ i^A = ^ '■ ^2 : —1 : 2, 

0, Qe = C2(5ci + 10c2 + 15c3 + 6C4), 

(48) 



1, Qi = 0, Qe = C2(5ci + 10c2 + 15c3 + 6C4), 

2, Qi =0, Qg = ci, (54 = (3ci + 2c4)pi, 

3, Q4 = 0, Qe = C4(5c3 + IOC4 + 15ci + 6C2), 

4, Oi - 0, 03== C3, Q4 == (3c3 + 2c2)pi. 



Since the coefficient pi must remain arbitrary, this leads to the Hst (21). 
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